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Abstract 



We compute the two-loop renormalization functions, in the RI' scheme, 
of local bilinear quark operators ipTip, where T denotes the Scalar and Pseu- 
doscalar Dirac matrices, in the lattice formulation of QCD. We consider both 
the flavor non-singlet and singlet operators; the latter, in the scalar case, leads 
directly to the two- loop fermion mass renormalization, Z m . 

As a prerequisite for the above, we also compute the quark field renor- 
malization, Z^p, up to two loops. 

We use the clover action for fermions and the Wilson action for gluons. 
Our results are given as a polynomial in csw, m terms of both the renormal- 
ized and bare coupling constant, in the renormalized Feynman gauge. We 
also confirm the 1-loop renormalization functions, for generic gauge. 

Finally, we present our results in the MS scheme, for easier comparison 
with calculations in the continuum. 

The corresponding results, for fermions in an arbitrary representation, are 
included in an Appendix. 

Keywords: Lattice QCD, Lattice perturbation theory, Fermion bilinears, 
clover action. 
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I. INTRODUCTION 



Studies of hadronic properties using the lattice formulation of QCD rely on the com- 
putation of matrix elements and correlation functions of composite operators, made out of 
quark fields. A whole variety of such operators has been considered and studied in numeri- 
cal simulations, including local and extended bilinears, and four-fermi operators. A proper 
renormalization of these operators is most often indispensable for the extraction of physical 
results from the lattice. 

In this work we study the renormalization of fermion bilinears O = ipTip on the lat- 
tice, where r = 1, 75. We consider both flavor singlet and nonsinglet operators. The cases 
T = 7/i, 75 7/i, 75 cr fJjU , will be presented in a sequel to this work. In order to obtain the renor- 
malization functions of fermion bilinears we also compute the quark field renormalization, 
Z^p, as a prerequisite. 

We employ the standard Wilson action for gluons and clover-improved Wilson fermions. 
The number of quark flavors Nf, the number of colors N c and the clover coefficient Csw are 
kept as free parameters. 

Our two-loop calculations have been performed in the bare and in the renormalized 
Feynman gauge. For 1-loop quantities, the gauge parameter is allowed to take arbitrary 
values. 

The main results presented in this work are the following 2-loop bare Green's functions 
(amputated, one-particle irreducible (1PI)): 

• Fermion self-energy: E^(g, a L ) 

• 2-pt function of the scalar ipip : T,g(qa L ) 

• 2-pt function of the pseudoscalar ip^ip : £p(ga L ) 

(a L : lattice spacing, q : external momentum) 

In general, one can use bare Green's functions to construct Z@' Y , the renormalization 
function for operator O, computed within a regularization X and renormalized in a scheme 
Y. 

We employ two widely used schemes to compute the various 2-loop renormalization 
functions: 

• The RI' scheme: Z^' RI , Z$ ,RI , Zp RI 

• The MS scheme: Z^ , Z^ , Z L P ^ 

The flavor singlet scalar renormalization function is equal to the fermion mass multi- 
plicative renormalization, Z m , which is an essential ingredient in computing quark masses. 

For convenience, the results for Z Q ' are given in terms of both the bare coupling constant 
g Q and the renormalized one: g^v , 9ms ■ 

Finally, as one of several checks on our results, we construct the 2-loop renormalized 
Green's functions in RI': S^ 7 (q, p,) (O = ip,S,P), as well as their counterparts in MS: 
T,Q S (q,fi). The values of all these functions, computed on the lattice, coincide with values 
computed in dimensional regularization (as can be inferred, e.g., from [1]). 
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The present work is the first two-loop computation of the renormalization of fermion 
bilinears on the lattice. One-loop computations of the same quantities exist for quite some 
time now (see, e.g., [2], [3], [4] and references therein). There have been made several 
attempts to estimate Z non-perturbatively; recent results can be found in Refs. [5-10]. 
Some results have also been obtained using stochastic perturbation theory [11]. A related 
computation, regarding the fermion mass renormalization Z m with staggered fermions can 
be found in [12]. 

The paper is organized as follows: Section II provides a formulation of the problem, along 
with all necessary definitions of renormalization schemes and of the quantities to compute. 
Section III describes our computational methods and the results which are obtained. Finally, 
in Section IV we discuss some salient features of our calculation, and comment on future 
extensions to the present work. 

Recently, there has been some interest in gauge theories with fermions in representations 
other than the fundamental. Such theories are being studied in various contexts [13-18], 
e.g., supersymmetry, phase transitions, and the 'AdS/QCD' correspondence. It is relatively 
straightforward to generalize our results to an arbitrary representation; this is presented in 
the Appendix. 

II. FORMULATION OF THE PROBLEM 
A. Lattice action 

We will make use of the Wilson formulation of the QCD action on the lattice, with the 
addition of the clover (SW) [19] term for fermions. In standard notation, it reads: 

Sl = S G + J2J2( 4r + m )j} f (x)i/} f (x) 
f x 



9 

* f x,p 



1 _ ^ 

+ 7 Csw 1>f( x ) °> F ^( x ) ^7( x )> (!) 



where : = — (Q^ - Q vlt ) (2) 

and . Q fiu U x , x+fi U x +fi, x+^+v Ux+^+u, x+u U x -\-v, x 

£^X, X+f Ux+V,X+V—IX U x -^i/ — jU; x — n U x — ^ x 
Uftx — n U x — n t x — [i — v U x — fj, — V, x — v U x — Vx 

~\~ U Xj x — v U x — v, x — u+fj, U x — v+n, x +n U x -\-n, x (3) 

Sg is the standard pure gluon action, made out of 1 x 1 plaquettes. The clover coefficient 
csw is treated here as a free parameter; r is the Wilson parameter (set to r = 1 henceforth); 
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/ is a flavor index; a^ v = (i/2)[y ll , 7^]. Powers of the lattice spacing a L have been omitted 
and may be directly reinserted by dimensional counting. 

The "Lagrangian mass" m Q is a free parameter here. However, since we will be using 
mass independent renormalization schemes, all renormalization functions which we will be 
calculating, must be evaluated at vanishing renormalized mass, that is, when m is set equal 
to the critical value m cr : m — > m cr = + O(gl). 



B. Definition of renormalized quantities 

As a prerequisite to our programme, we will need the renormalization functions for the 
gluon, ghost and fermion fields (A®, c a , ip), and for the coupling constant g and gauge 
parameter a, defined as follows: 




g Q = // Z g g, a Q = Z a 1 Z A a (4) 

The value of each Zq depends both on the regularization X and on the renormalization 
scheme Y employed, and thus should properly be denoted as Z G ' . The scale \i enters the 
relation between g Q and g only in dimensional regularization (D = 4 — 2e dimensions). 

We will need Za, Z c , Z a wAZ g to 1 loop and Z^ to 2 loops. Our 1-loop results, per- 
formed in a generic gauge, are in agreement with results found in the literature (see, e.g., 
Refs. [4,20]). 



C. Definition of the RI' scheme 



This renormalization scheme [21-23] is more immediate for a lattice regularized theory. 
It is defined by imposing a set of normalization conditions on matrix elements at a scale /2, 
where (just as in the MS scheme) [24]: 

At = /i(47r/e^) 1 / 2 (5) 

(7e is the Euler constant). 

In Euclidean space, the fermion self energy S^(g, a L ) = ii + m Q + O(g^) is renormalized 
through: 



lim 

a L ^0 



4^'(a L /i)tr(s^g,a L )^/(^g 2 



(6) 



The trace here is over Dirac indices; a Kronecker delta in color and in flavor indices has 
been factored out of the definition of 



Similarly, for the ghost self energy S^(g, a L ) —q+ 0{g ): 



lim 



'c 1°lMJ 75 



(7) 
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Za and Z a are extracted from the gluon propagator G^ u (q,a L ) with radiative 



corrections 1 : 



n T (a L g) 



n L (a L g) 



(8) 



where n TjL (a L g) = 1 + O(gl). The normalization conditions are: 



lim 



lim 



tL,RI' 



(a L /i) 



n T (a L g) 
1 

n L (a L g) 



= 1 



q 2 =fi 2 



(9) 
(10) 



We have checked explicitly that Z^' RI ' = 1 up to one loop, in agreement with the continuum. 

For consistency with the Slavnov- Taylor identities, Z g in the RI' scheme is defined as in 
the MS scheme. In dimensional regularization {DR) this is achieved by requiring that the 
gluon-fermion-antifermion 1PI vertex function, G A ^, renormalizes as follows [1]: 



lim 



/'"-(finite 



(11) 



The value of Z g is tuned in such a way as to absorb only the poles in e which appear in 
G A ^ (together with matching powers of ln(47r) — 7e); this leads to a result for G A n ^ which 
is finite but not unity. Before rescaling, we have first divided G A ^ by the bare coupling 
constant, as in Ref. [25], in order to have unity as the tree level value for G A n ^. We have set 
the fermion momentum to zero; q refers to the gluon/ ant if ermion momentum. Alternatively, 
a similar procedure can be performed on the gluon-ghost-antighost vertex: 



lim 



z dr,rv {z n R , R i'y/2 Z n R ,Ri' GAM 



/'-(finite 
^Acc 



(12) 



Eq.(12) leads to exactly the same value for Z g . 

The corresponding renormalization conditions on the lattice read: 



lim 

a L ^0 



Zi' RI '(Z^ 2 Z^'G^( q ,aJ 



g 2 =fi 2 



or, equivalently: 



lim 



Zt RI '(ZA> R y /2 Zi> RI 'G L Mc (q,a L 



/'-(finite 

u 'Ai/n/j 



/'-(finite 
^Acc 



(13) 



(14) 



where the expressions G A n ^ and G A n £ e are required to be the same as those stemming from 
the continuum (Eqs.(ll),(12)). We have calculated Z^ ,RI \ using either one of Eqs.(13), (14), 
and have verified that the same result is obtained. 



One should carefully distinguish among the following standard symbols: a L : lattice spacing; a a , 
a Ri' i a M~9 : bare and renormalized gauge parameters. 
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D. Conversion to the MS scheme 



For easier comparison with calculations coming from the continuum, we need to express 

L RI ! 

our results in the MS scheme. Each renormalization function on the lattice, Z ' , may be 
expressed as a power series in the renormalized coupling constant gRi>. For the purposes of 
our work the conversion of gRp to MS is trivial since: 

9Ri' = 9ms + °(9ws) (15) 

As already mentioned, our 1-loop calculations for Z^, Z c , Z a and Z g are performed in a 
generic gauge, cxrp. The conversion to the MS scheme is given by [26]: 

a Ri' = V L,RP a MS = "Is/ Ca(a^ms, «1s) ( 16 ) 

Since the ratio of Z's appearing in Eq.(16) must be regularization independent, it may be 
calculated more easily in dimensional regularization [1]; to 1 loop, the conversion factor Ca 
equals: 

r,DR,RI' 2 

C A (g, a) = = 1 + ^f^y [(9a 2 + 18a + 97) N c - 40N f ] (17) 



J A 



(Here, and throughout the rest of this work, both g and a are in the MS scheme, unless 
specified otherwise.) 

Once we have computed the renormalization functions in the RV scheme we can construct 
their MS counterparts using conversion factors which, up to the required perturbative order, 
are given by: 
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yh,RV yDR,RI' 

CJg, a) = -£= = c _ = 1 + N c (18) 

cVy ' ' ^L,MS ^DR,MS 1671-2 c v 1 



C^{g,a) 



rvl^^RI 1 ryTDR^Rl 1 



7 L,MS 7 DR,MS 

g 2 g A 

1 o Cp a -\ ; — - Cp 

16tt 2 * 8(16tt 2 ) 2 



(8a 2 + 5) c F + 14 Nf 



(9a 2 - 24C(3) a + 52a - 24((3) + 82) N c 



(19) 



where cp = (iV 2 — 1)/(2N C ) is the quadratic Casimir operator in the fundamental rep- 
resentation of the color group; ((x) is Riemann's zeta function. (We employ a standard 
normalization for the generators of the algebra, T a , see the Appendix.) 

E. Renormalization of fermion bilinears 

The lattice operators Or = ipTif) must, in general, be renormalized in order to have 
finite matrix elements. We define renormalized operators by 
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Of = Z^ RI '(a L fi)O ro 



(20) 



The flavor singlet scalar operator receives also an additive renormalization, which must 
be taken into account; we discuss this issue in the following subsection. For the scalar 

L RI 1 

(S) and pseudoscalar (P) operators, the renormalization functions Z T ' can be obtained 
through the corresponding bare 2-point functions Ef (qa L ) (amputated, 1PI) on the lattice, 
in the following way: 



lim 

a L ->0 

lim 



7L,RI' Z L,RI' E L {qa j 



yL,RI ! rvl^^RV 



E£(ga L ) 



1 

75 



(21) 
(22) 



where: 



E§(ga L ) = U + C% 2 ) , E£(ga L ) = 75 + 0(g 2 ) (23) 

L RI 1 

Once the quantities Z T ' have been calculated, one may proceed to compute them also 
in the MS scheme. In the case of the scalar operator (Os Q = 4>o4>o), the renormalization 



function, Zg' MS ', can be obtained by: 



rL,MS 



Z L S > RI '/C s (g,a) 



(24) 



where Cs(g, ol) is a regularization independent conversion factor and has been calculated in 
dimensional regularization [1]: 



i^L^RI 1 



7 DR,RV 



V DR,MS 
Z 5 



1+ lfe CF (a + 4) + 
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cf 



(24a 2 + 96a - 288C(3) + 57) c F + 166 N f 

(25) 



24(16vr 2 ) 2 

- (l8a 2 + 84a - 432((3) + 1285) N t 



The treatment of the pseudoscalar operator (Op Q = V^TsVO m the MS scheme requires 
special attention, due to the non-unique generalization of 75 to D dimensions. A practical 
definition of 75 for multiloop calculations, which is most commonly employed in dimensional 
regularization and does not suffer from inconsistencies is [27]: 



. 1 

75 — % vi 1/3 U4 7^1 7^2 7^3 7^4 



v l = 0, 1, 2, 3 



(26) 



Of course, 75 as defined in Eq.(26) does not anticommute with the D-dimensional 7^; an 
ultimate consequence of this fact is that Ward identities involving the axial and pseudoscalar 
operators, renormalized in this way, are violated. 

To obtain the correctly renormalized pseudoscalar operator, one must introduce an ex- 
tra finite factor, Z 5 , in addition to the usual renormalization function Zp R,MS which only 
contains poles in e. We set: 
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P = Z 5 (g)Z° R > MS Po (27) 

Z 5 is defined by the requirement that the scalar and pseudoscalar renormalized Green's 
functions coincide: 



(2g) 

Gp S 

The value of Z 5 , calculated in Ref. [28], is gauge independent; it equals: 

z ^ = 1 - 1& (8 CF> + j^fih N ' + h N >) + (29) 

Zp MS can now be obtained by: 

Z%& 3 = Z£ m '/(C S Z 5 ) (30) 

Similarly, one can convert the RI' renormalized Green's functions, G^ 1 ' , to their MS 
counterparts, through: 

qRV qRI' 

7m = c * Cs ' rh = c ^ CsZ5 ( 31 ) 

(In Eqs.(30, 31) it is understood that powers of gni>, ctRv, implicit in RI' quantities, must 
also be converted to gjjg, otj^g, respectively, using Eqs.(15, 16). ) 



F. Fermion Mass Renormalization 

As a by-product of this work, one can evaluate the fermion multiplicative mass renor- 
malization, Z m , which is directly related to the scalar flavor singlet operator. This operator 
differs from the ones considered thus far, in that it receives also an additive renormalization, 
since it has a nonzero perturbative vacuum expectation value; thus, it mixes with the iden- 
tity at the quantum level. Once its vacuum expectation value is subtracted, the resulting 
operator is multiplicative renormalizable. The renormalization is then simply given by only 
connected diagrams of the original operator (Figs. 3, 4 and 5); all disconnected diagrams 
are easily shown to cancel out. 

The perturbative vacuum expectation value 2 is of course a power divergent quantity, and 
it cannot be expected to approach well the value of the corresponding disconnected matrix 
elements in numerical simulations. Fortunately, this quantity is not needed for multiplicative 
renormalization, as mentioned above. However, as regards simulations, one should bear in 
mind that disconnected parts must be evaluated and subtracted from matrix elements, before 
the latter can be renormalized. 



2 For a tree level computation of this quantity, see Ref. [29]. 
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Let us express the fermion self energy in the following way: 



^lb RI = H^odd (ga L ,mod L ,3o) + 1 • — S cven (ga L ,m a L ,5f ) 



(32) 



o a L + C^o)- Terms like <lfi 7m/ though a priori 



where E odd = 1 + 0{gl) and S cvcn = m, 
allowed by hypercubic symmetry, are eventually seen to cancel, as expected by Lorentz 
invariance. 

For generic values of m , the even part of E^' is power divergent; in order to achieve a 
finite renormalized mass, m r , the values of the Lagrangian mass m must be near a critical 
value, m cr , at which E cvcn vanishes: T* even (qa h , m cr a L , g Q ) = + 0(q 2 a 2 ). That is, m cr is 
required to satisfy: 



E eV en(0, m cr a L , g Q ) = 



(33) 



This is a recursive equation which can be solved for m cr order-by-order in perturbation 
theory. Its value is known to two loops for Wilson fermions: [30] (confirmed independently 
in [31]), and for clover fermions: [32] (with Wilson gluons), [33] (with Symanzik gluons). 
Only the 1-loop value of m cr enters the present calculation. 

We can perform a Taylor expansion with respect to the bare mass 3 , m B = m Q — m CT , for 
both E odd and E even : 



S dd(ga L ,m a L ,5f ) = 



S dd(ga L ,m a L ,5( ) 



m =m CT 



+ 0{m B a L ) 



(34) 



" E eve n W a L > m o<\, <? ) 



+ m B 



Eeven (?° L ' m o a L J 9o 

d 



d(m a L ) 



^cvcn 



+ 0(a L 



(35) 



to — i/tcr 



Note that when m = m cr , the first term on the r.h.s. of Eq.(35) vanishes in the limit 
a L — > 0, by virtue of Eq.(33). 



Having in mind that, in calculating E , one is interested in the limit a L — > 0, the 



fermion self energy takes the form: 



+ 1 • m B 



E dd(ga L ,m a L ,5( ) 
9 



to — "<<cr 



_<9(m aJ 

The renormalized fermion mass is now defined by: 



^cvcn 



(36) 



to— "tcr 



3 Note that m cr (and, consequently, m ) is power divergent in a L since its calculation contains no 
other dimensional quantities; m B , on the other hand is at most logarithmically divergent in a L . 
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The renormalization condition for Z^ ,RI (Eq.(6)) for nonzero m r , reads: 



lim 



[id + m r 



= 



(37) 



(38) 



By combining Eqs.(37) and (38) we find the renormalization condition for Z, 

d r 



L,RI'. 
m 



lim 



7 L,RI' yL.RI' 



d(m a h )' 



,(ga L ,m a L ,^ ) 



m =m CT 



(39) 



We stress again that, even though the Lagrangian mass m may take arbitrary values, the 
renormalization condition involves only m — > m CT . 

In order to establish a relation between Z^ RI ' and Z^f^ let , note that Eq.(39) coincides 
with Eq.(21) if d/d (m a L )E even = ^s'^Lgiet- Indeed, the equality between d/d (m a L )E even 
and ^sl'singiet holds diagram by diagram in perturbation theory, noting that: 

• The tree level value equals 1, in both cases 

• The effect of inserting the scalar operator on a given fermion propagator of any 
self-energy Feynman diagram is equivalent to taking the negative partial derivative 
—d/d(m a L ) of that propagator 

• Combinatorial factors agree 

• There is an extra minus sign in the geometric series summation of 1PI diagrams leading 
to the fermion self-energy 

Once all of the above statements are taken into account, one comes to the conclusion 
that: 



yL,RV _ yL,RV 
m S, singlet 



(40) 



Given that m cr a L = 0(g%), all two-loop calculations can be performed with strictly 
massless fermion propagators, provided that appropriate fermion mass counterterms are 
introduced on one-loop diagrams. 



III. COMPUTATION AND RESULTS 

The Feynman diagrams relevant to the fermion self-energy £^(g, a L ), at 1- and 2-loop 
level, are shown in Figs.l and 2, respectively; those relevant to S^(ga L ), Sp(ga L ) are shown 
in Figs. 3 and 4. 

For flavor singlet bilinears, there are 4 extra diagrams, in addition to those of Fig.4, 
shown in Fig.5; in these diagrams, the operator insertion occurs inside a closed fermion 
loop. 
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The evaluation and algebraic manipulation of Feynman diagrams, leading to a code 
for numerical loop integration, is performed automatically using our software for Lattice 
Perturbation Theory, written in Mathematica. 

The most laborious aspect of the procedure is the extraction of the dependence on the 
external momentum q. This is a delicate task at two loops; for this purpose, we cast algebraic 
expressions (typically involving thousands of summands) into terms which can be naively 
Taylor expanded in q to the required order, plus a smaller set of terms containing superficial 
divergences and/or sub divergences. The latter can be evaluated by an extension of the 
method of Ref. [34] to 2 loops; this entails analytical continuation to D > 4 dimensions, 
and splitting each expression into a UV-finite part (which can thus be calculated in the 
continuum, using the methods of Ref. [35]), and a part which is polynomial in q. A primitive 
set of divergent lattice integrals involving gluon propagators, which can be obtained in this 
manner, can be found in Ref. [36]. 

Some of the diagrams contributing to S^(g, a L ), S^(ga L ) and T,p(qa L ) are infrared diver- 
gent when considered separately, and thus must be grouped together in order to give finite 
results. Such groups are formed by diagrams (7-11), (12-13), (14-18), (19-20), (21-23) in 
Fig. 2, diagrams (3-7), (8-9), (10-11,19) in Fig. 4 and diagrams (1-2), (3-4) in Fig. 5. 

In Figures 1-5, "mirror" diagrams (those in which the direction of the external fermion 
line is reversed) should also be taken into account. In most cases, these coincide trivially 
with the original diagrams; even in the remaining cases, they can be seen to give equal 
contribution, by invariance under charge conjugation. 

As mentioned before, all calculations should be performed at vanishing renormalized 
mass; this can be achieved by working with massless fermion propagators, provided an 
appropriate fermion mass counterterm is introduced (diagram 23 in Fig. 2 and diagram 11 
in Fig.4). 

All two-loop diagrams have been calculated in the bare Feynman gauge (a Q = 1). One- 
loop diagrams have been calculated for generic values of a D ; this allows us to convert our 
two-loop results to the renormalized Feynman gauge (am' = 1 or aj^g = 1). 

Numerical loop integration was carried out by our "integrator" program, a metacode 
written in Mathematica, for converting lengthy integrands into efficient Fortran code. Two- 
loop numerical integrals were evaluated for lattices of size up to L = 40; the results were 
then extrapolated to L — > oo. Extrapolation is the only source of systematic error; this error 
can be estimated quite accurately (see, e.g. Ref. [37]), given that L-dependence of results 
can only span a restricted set of functional forms. 
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FIG. 2. Two-loop diagrams contributing to Z^. Wavy (solid, dotted) lines represent gluons 
(fermions, ghosts). Solid boxes denote vertices stemming from the measure part of the action; a 

solid circle is a fermion mass counterterm. 
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FIG. 3. One-loop diagram contributing to Z$ and Zp. A wavy (solid) line represents gluons 
(fermions). A cross denotes the Dirac matrices II (scalar) and 75 (pseudoscalar). 




1 2 3 4 5 6 




7 8 9 10 11 12 




13 14 15 16 17 18 




19 20 

FIG. 4. Two-loop diagrams contributing to Z$ and Zp. Wavy (solid, dotted) lines represent 
gluons (fermions, ghosts). A solid box denotes a vertex from the measure part of the action; a 
solid circle is a mass counterterm; crosses denote the matrices II (scalar) and 75 (pseudoscalar). 




12 3 4 

FIG. 5. Extra two-loop diagrams contributing to Z$, singlet- A cross denotes an insertion of a 
flavor singlet operator. Wavy (solid) lines represent gluons (fermions). 
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A. One-loop results 



1-loop results for Z^' RI , Zg' RI and Zp' RI are presented below in a generic gauge. The 
errors result from the L — > oo extrapolation. 



rvL.RI 1 -i . 9o 



n 2 

Zq = 1 + n Cp 



16tt 2 



n 2 

rjL,RV -i . i/o 



167T 2 



(ln(aj/i 2 ) - 4.792009570(1)) a G + 16.644413858(5) 
-2.248868528(3) c sw - 1.397267102(5) c 2 sw 

3 ln(aj/x 2 ) - a Q - 16.9524103(1) 
-7.7379159(3) c sw + 1.38038065(4) c§ w 

3 ln(a 2 /i 2 ) ~ «o - 26.5954414(1) 
+2.248868528(3) c sw - 2.03601561(4) c 2 sw 



The corresponding quantities in the MS scheme are: 

(ln(aj/i 2 ) - 3.792009570(1)) a + 16.644413858(5) 



zi*™ = i + jL CF 



16tt 2 



-2.248868528(3) c sw - 1.397267102(5) c 



sw 



Z L S MS = 1 + 



,L,MS 



= 1 + 



16tt 2 



16tt 2 



c f 



3 ln(a 2 /z 2 ) - 12.9524103(1) 

-7.7379159(3) c sw + 1.38038065(4) c| w 
3 ln(a 2 /i 2 ) - 14.5954414(1) 

+2.248868528(3) c sw - 2.03601561(4) c| w 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



Our results confirm the existing results found in the literature [4] (note, however, a 
difference in Zp MS ; this is entirely due to the factor Z 5 in Eq.(30) ). 
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B. Two-loop results 



The evaluation of all Feynman diagrams in Figs. 1-5 leads directly to the corresponding 
bare Green's functions E^, E^ and Ep. These, in turn, can be converted to the corresponding 



renormalization functions Z^ Y , and Zp' Y (Y = RI' or MS), via Eqs.(6), (21) and (22). 
To this end, we need the following one-loop expression for Z^ (recall that Z a — lto this 
order) : 



r L,RI' 



,L,MS 



= 1 + 



16tt 2 



In (a 2 fx 2 



-N f --N c 
3 i 3 



+N f (-2.168501047(1) + 0.7969452308(4) c sw - 4.7126914428(1) c 



sw 



+39.47841760436(1) c F + 1.94017130069(1) N c 



(47) 



To express our results in terms of the renormalized coupling constant, we also need the 



one-loop expression for Z g 



r?L,RI' _ r?L,MS 



16tt 2 



+ O (gt) 
In (a 2 ft 



--N f + — N c 
3 } 6 



+A^/ (0.5286949677(5) - 0.3984726154(2) c sw + 2.35634572140(7) d? w ) 



•19.73920880218(1) c F - 3.54958342046(1) N c 



+ O (gi) 



(48) 



Eqs.(47, 48) are in agreement with older references (see, e.g., Ref. [20]). 

L RI 1 L RI 1 L RI 1 

We present below Z^ , Z s ' and Zp to two loops in the renormalized Feynman 

gauge am' = 1; we also present the MS analogues Z^' MS , Z$' MS and Zp MS in the gauge 
= 1. For conciseness, we omit the bare Green's functions; it is a straightforward exercise 
to recover these from the corresponding Z's. 
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7 L,RP 



1 + 



+ 



16tt 2 
9t 



c F 



\n(a 2 h 2 ) + 11.852404288(5) - 2.248868528(3) c sw - 1.397267102(5) c 2 w 



c F 



(16 ^ L ta W)(^+^/-^) 

+ ln(aj/i 2 ) ( - 6.36317446(8) N f + 0.79694523(2) N f c sw 
-4.712691443(4) N f c| w 
+49.83082185(5) c F - 2.24886861(7) c F c sw 
-1.39726705(1) c F c 2 sw + 29.03029398(4) iV c ) 
+N f ( - 7.838(2) + 1.153(1) c sw + 3.202(3) c 2 sw 

+6.2477(6) 4 W + 4.0232(6) 4 W ) 
+c F (505.39(1) - 58.210(9) c sw + 20.405(5) c 2 sw 

+18.8431(8) c 3 sw + 4.2793(2) c* w ) 
+N C ( - 20.59(1) - 3.190(5) c sw - 23.107(6) c 2 sw 

-5.7234(5) c 3 sw - 0.7938(1) c 4 sw ) 



(49) 



Z L,MS = 1 + 



16tt 2 

gt 



C F 



\n(a 2 L jj, 2 ) + 12.852404288(5) - 2.248868528(3) c sw - 1.397267102(5) 4 W 



c F 



ln 2 (a 2 /, 2 ) (lNf + lc F - 8 -N c 



(16tt 2 ) 2 

+ ln(a 2 /i 2 ) ( - 4.58539668(8) JV/ + 0.79694523(2) N f c sw 
-4.712691443(4) N f 4 w 
+50.83082185(5) c F - 2.24886861(7) c F c sw 
-1.39726705(1) c F 4w + 21.91918287(4) Nc) 
+N f ( - 15.970(2) + 1.950(1) c sw - 1.510(3) c 2 w 

+6.2477(6) 4 W + 4.0232(6) c| w ) 
+c F (556.10(1) - 60.459(9) c sw + 19.007(5) c§ w 

+18.8431(8) 4 W + 4.2793(2) c^ w ) 
+N C (13.68(1) - 3.190(5) c sw - 23.107(6) 4 

-5.7234(5) 4 W - 0.7938(1) c| w ) 



.2 



(50) 



16 



1 + 



+ 



16tt 2 
9t 



c F 3 \n(a 2 fi 2 ) - 17.9524103(1) - 7.7379159(3) c sw + 1.38038065(4) c 2 sw 

(l 67 rT CF [ lnW) (l CF + Nf -^ Nc ) 

+ ln(a^ 2 ) ( - 8.1721694(5) N f + 2.3908354(3) N f c sw 



■14.13807433(4) N f q 



sw 



+66.0780218(9) c F - 23.213749(2) c F c sw 
+4.1411425(3) c F c 2 sw + 55.7975008(9) iV c ) 
+N f (24.003(3) + 11.878(5) c sw + 25.59(1) c 2 sw 

+22.078(3) c 3 sw - 6.1807(8) c| w ) 
+c F ( - 602.35(6) + 91.07(7) c sw + 51.15(5) c| 

-27.759(4) c 3 sw - 2.688(1) c^ w ) 
+iV c ( - 38.16(4) - 132.40(5) c sw - 4.04(3) c§ 



.2 

-sw 



.2 

-sw 



+ 12.576(3) ci w + 1.0175(8) 4 W 



(51) 



1 j 9o 

1 + TTT^T C F 



+ 



167T 2 



'sw 



(16tt 2 ) 



■c F 



3 ln(aj/i 2 ) - 12.9524103(1) - 7.7379159(3) c sw + 1.38038065(4) 4 
ln 2 (a 2 ^ 2 ) (iy / + ^c F -^iV c ) 

+ ln(a 2 /i 2 ) ( - 4.8388361(5) TV/ + 2.3908354(3) N f c sw 
-14.13807433(4) N f 4 w 
+81.0780218(9) c F - 23.213749(2) c F c sw 
+4.1411425(3) c F e| w + 37.4641674(9) iV c ) 
+N f (10.688(3) + 15.863(5) c sw + 2.02(1) 4w 

+22.078(3) 4 W - 6.1807(8) c| w ) 
+c F ( - 462.67(6) + 52.38(7) c sw + 58.05(5) c| w 

-27.759(4) 4 W - 2.688(1) c| w ) 
+iV c (36.93(4) - 132.40(5) c sw - 4.04(3) c§ w 

+12.576(3) 4 W + 1.0175(8) 4 W ) (52) 
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1 + TTT^T C F 



+ 



16tt 2 

at 



3 ln(a 2 /T) - 27.5954414(1) + 2.248868528(3) c sw - 2.03601561(4) c 2 sw 

(16 , 2)2 ^[ln 2 (a^ 2 ) (lc F + N f -^N c ) 

+ ln(aj/i 2 ) ( - 8.1721694(4) N f + 2.39083540(6) N f c sw 



■14.13807433(4) N f c : 



sw 



+37.1489292(7) c F + 6.746606(1) c F c sw 
-6.1080465(3) c F c 2 sw + 55.7975008(7) JV C ) 

+N f (38.231(3) - 7.672(5) c sw + 55.32(1) c 2 sw 

-7.049(3) c 3 sw + 4.7469(8) c^ w ) 
+c F ( - 876.98(4) + 84.52(2) c sw + 38.65(4) c| w 

+19.974(3) c 3 sw + 2.873(1) c^ w ) 
+iV c ( - 104.35(3) - 38.06(2) c sw - 14.57(3) c 2 sw 

-4.429(2) 4 W - 1.2898(7) 4 W 



(53) 



1 J 0° 

1+ 16^ 

,4 



3 ln(a 2 /T) - 14.5954414(1) + 2.248868528(3) c sw - 2.03601561(4) c 2 sw 
c F 



2, 2 -2\ ! "\ " 9 „ ' ' y 



ln>^)(7V / + _ CF - 2 



(16tt 2 ) 2 

+ ln(a2/i 2 ) (0.4944972(4) TV/ + 2.39083540(6) iV> c sw 
-14.13807433(4) TV/ 
+76.1489292(7) c F + 6.746606(1) c F c sw 
-6.1080465(3) c F c| w + 8.1308341(7) N c ) 

+N f (16.013(3) + 2.688(5) c sw - 5.94(1) c 2 sw 

-7.049(3) 4 W + 4.7469(8) c| w ) 
+c F ( - 586.45(4) + 113.76(2) c sw + 12.18(4) c 2 sw 

+19.974(3) c 3 sw + 2.873(1) c| w ) 
+iV c (27.31(3) - 38.06(2) c sw - 14.57(3) 

-4.429(2) 4 W - 1.2898(7) 4 W ) 



-sw 



(54) 



18 



All expressions reported thus far for Z$ and Zp refer to flavor non singlet operators. In 
the case of Zp, all diagrams of Fig. 5 vanish, so that singlet and non singlet results coincide, 
just as in dimensional regularization. For Z s on the other hand, the above diagrams give 
an additional finite contribution: 



7 L,RI' 
J S, singlet 



7 L,RI' 



9t 



(16tt 2 ) ; 



cpNf ( - 107.76(1) + 82.27(2) c sw - 29.727(4) c; 



.2 

sw 



+ 3.4400(7) 4 W + 2.2758(4) 4 W 



(55) 



The same extra finite contribution applies also to the MS scheme. 

Finally, for completeness, and as an additional check on our results, we compute the 
renormalized Green's functions (for v anishing renormalized mass): 



G™'(fi/ q ) = Zi> RI 'x 



iRI' i 
iRI' 



odd 

L III' ryL,Rr yj. 



G f\n/ q ) = z^ z, 

Gfm/ q )^z^'z L P ^ 



yL 



(56) 
(57) 
(58) 



Similarly for MS, taking into account Eq.(31). 

Since these functions are regularization independent, they can be calculated also using, 
e.g., dimensional regularization. We have computed G^, Gs and Gp in both ways: either 
starting from our Eqs. (47-54) or using renormalization functions from dimensional regular- 
ization [1]. In all cases the two ways are in complete agreement. We obtain: 



+ 



16tt 2 
9rv 
(16tt 2 )- 



c F 



+ H?/ q ')[-^N f -\cp + 2 -^N c 



" ' f 16^ CF 

(16vr 2 )- 



+ 



Cp 



InV/V) 



c F + N c 



+ H?/ q ')[-N f - l -c F + 2 ±N c 



7 Ar 5 

— Nt c F 

4 } 8 



/143 



6C(3) N, 



(59) 



(60) 



1 _l_ g2RV r. 



4 ln(^) 



, 9rp 

+ ONO CP 



(16tt 2 ) ; 



InV/V) [~Nf + 8 c F + -| 7V C ) 



(61) 



Eq.(61) holds also for the case of the pseudoscalar operator: Gp 1 ' = Gg 1 ' . 
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Q— 



(16tt 2 )' 



4 ln(/i 2 /g 2 ) + 6 

ln 2 (/i 2 /g 2 ) (-N f + 8c F + ^-N t 
+ ln(/, 2 /g 2 ) (-^ iV, + 24 c F + ^ iV ( 
+ (-yiy /+ (22 + 12C(3))c F 



c 

227 



24C(3) N ( 



(62) 



qMS 



I _l_ 9ms 



16tt 2 

4 

MS 



CF 



(167T 2 ) 2 



4 ln(/i 2 /g 2 ) + 14 

InV/g 2 ) (-iV / + 8c F + ^iV c ) 



16 v , , cc „ , 130^ 



+ \nm 2 /q 2 )[--N f + 5Qc F + ., 



9 



iV> + ( 134 + 12C(3) ) CF +(™- 24C(3)) 7V C ) 



In Figs. (6a, 6b), (7a, 7b), (8a,8b) we plot 



t 7 L,MS 



V 9 

zi RI \ (z 



L,MS rvL,RV\ 
S i ^S i 



and (Z 



(63) 



L,MS 



7 L,MS 
J S > 



7 L,MS , 7 L,MS 
Zj p ana A s, singlet 



N f = 2 and /3 Q = 5.3. The 



Z P ' ), respectively, as a function of csw- 111 practice, of course, only specific values of csw 
are relevant, in the range 1 < csw < 1-8, corresponding to perturbative or non-perturbative 
determinations. For definiteness, we have set iV c = 3, Jx = l/a L and (3 Q = 2N c /g1 = 6.0. Our 
results up to two loops for each Z are shown for both Nf = and iVj = 2, and compared 
to the corresponding one-loop results. Furthermore, in the scalar case, we also present the 
two-loop result for the flavor singlet operator. 

In Fig. 9 we present, on the same plot, the values of Z^ ,MS 
up to 2 loops, versus csw- We have chosen N c — 3, p, — l/a L 
corresponding results in the RI' scheme are plotted in Fig. 10. 

There are a number of non-perturbative (NP) estimates of renormalization constants 
in the literature, in the RI' scheme (see, e.g., [38-40]) and in the Schrodinger functional 
scheme [41]. Our 2-loop results still differ from NP results in RI', and this leaves open the 
possibility that higher loop effects may still be important, even though the perturbative 
series shows reasonable signs of convergence. A putative reason for this difference is the 
fact that the bare coupling constant g is known not to be a good expansion parameter. 
One may also express the renormalization functions in terms of the renormalized couplings: 
9jTs or 9R 1 '- T ne resu lti n g expressions for Z^ ,Y , Zp Y and Zg ,Y inglet (Y = MS , RI') as a 
function of csw> are shown in Figs. 11 and 12, for the same values for N c , Nf, p, and j3 as in 
Figs. 9 and 10. For values of the clover parameter beyond its typical range, csw > 1-8, the 
behavior of the renormalization functions shows signs of instability at the scale ft = l/a L . 
There exist also several alternative definitions of an effective coupling in the literature; one 
should be aware, however, that the use of many of these definitions (coming, e.g., from 
boosted perturbation theory) can only be justified for 1-loop quantities, not beyond. For 
this reason, we have preferred to provide the bare results in this paper, leaving to the reader 
the straightforward task of converting these results to their favorite scheme. 
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FIG. 6b. zh ,RI (a L /u) versus csw (N c = 3, \i = l/a L , (3 Q = 6.0). Same notation as in FIG. 6a. 
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c sw 



1.5 
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FIG. 7a. Zg' M (o L At) versus csw (N c = 3, p, = l/a L , /? = 6.0). Results up to 2 loops, for the 
flavor non-singlet operator, are shown for Nf = (dashed line) and Nf = 2 (solid line); 2-loop 
results for the flavor singlet operator, for Nf = 2, are plotted with a dash-dotted line; one-loop 

results are plotted with a dotted line. 
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FIG. 8a. Zp' MS (a L p,) versus csw (N c = 3, \i = l/a L , /3 = 6.0). Results up to 2 loops are shown 
for iVj = (dashed line) and Nf = 2 (solid line); one-loop results are plotted with a dotted line. 
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FIG. 8b. Zp' RI (a L/ u) versus csw (N c = 3, \i = l/a L , /3 Q = 6.0). Same notation as in FIG. 
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FIG. 9. Z^ ,MS (dotted line), Zg' M (dashed line), Zp' M (solid line) and Zg'^ nglet (dash-dotted 
line) up to 2 loops, versus csw (N c = 3, p = l/a L , Nf = 2, (5 Q = 5.3). 
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FIG. 10. Z^' RI (dotted line), Zg' RI (dashed line), Zp RI (solid line) and Zg ,R f nglet (dash-dotted 
line) up to 2 loops, versus csw (N c = 3, \x = l/a L , iVf = 2, /3 Q = 5.3). 
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FIG. 11. 2-loop renormalization functions Z s ' (dashed line), Zp' MS (solid line) and singlet 
(dash-dotted line) expressed in terms of the renormalized coupling constant (fjjg-, versus csw 



(N c = 3, fit = l/a L , N f = 2,f3 = 5.3). 
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FIG. 12. 2-loop renormalization functions Zg' RI (dashed line), Zp' RI (solid line) and Zg'^f let 
(dash-dotted line) expressed in terms of the renormalized coupling constant gui 1 , versus csw 

(N c = 3,p = l/a h1 Nf = 2 1 P = 5.3). 
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IV. DISCUSSION 



As can be seen from Figs.6a-8b, all 2-loop renormalization functions differ from 1-loop 
values in a significant way; this difference should be taken into account in MC simulations, 
in order to reduce systematic error. At the same time, 2-loop contributions are consistently 
smaller than 1-loop contributions, indicating that the (asymptotic) perturbative series are 
under control. 

The dependence on the clover parameter csw is also quite pronounced. In the present 
work, csw was left as a f ree parameter; its optimal value, as dictated by 0(a L ) improvement, 
has been estimated both non-perturbatively [42] and perturbatively (to 1-loop) [19]. 

Our results regard both the flavor nonsinglet and singlet operators. For the pseudoscalar 
operator, these cases coincide, just as in dimensional regularization. The scalar operator, 
on the other hand, receives an additional finite (a L /x independent) contribution in the flavor 
singlet case. Z s> S i ng i e t is seen to be equal to the fermion mass renormalization Z m , which is 
an essential ingredient in the computation of quark masses. 

We note also that, in dimensional regularization, both the scalar and pseudoscalar flavor 
singlet operators renormalize in the same way as their non-singlet counterparts, for mass 
independent renormalization schemes. Consequently, the conversion factors Cs and Cp, as 
well as Z 5 , stay the same for flavor singlets. 

A breakdown of our results on a per diagram basis has not been presented here, due to 
lack of space; it is available from the authors upon request. 

The 2-loop computation of the renormalization functions for the Vector, Axial and Tensor 
bilinears is work currently in progress. 

Besides the strictly local definitions of fermion bilinears, ipTip, one can consider a family 
of more extended operators (see, e.g., [42]), with the same classical continuum limit, as 
dictated by 0(a L ) improvement. The renormalization of these extended operators involves 
more Feynman diagrams, since their vertices may also contain gluon lines; however, the 
computation is actually less cumbersome, since all additional contributions are now free of 
superficial divergences. We will be reporting the results of this computation in a future 
work. 



26 



APPENDIX: Fermions in an arbitrary representation 

Our results for Zs, Zp, Eqs.(49, 51, 53), can be easily generalized to an action with 
Wilson/clover fermions in an arbitrary representation R, of dimensionality dp . 

In this case, the gluon part of the action remains the same, while all link variables 
appearing in the fermion part of the action assume the form: 

U x>x+l , = ex P (ig A«(x)T a ) — > U x , x+Ii = exp(i g Afa) T a R ) (64) 

Using standard notation and conventions, the generators T a in the fundamental representa- 
tion satisfy: 

/V 2 — 1 1 
[T a , T b ] = i f abc T c , T " T<1 = 11 °f = 1 c , tr(T a T b ) = 5 ab t F = 5 ah - (65) 

a 2iV c 2 

In the representation R we have: 

[T°,T b R ]=if abc T c R , Y. t rTp = Kc r , U{T a R T b R ) = 5 ab t R (66) 

a 

where: t R = {dp c R )/(N 2 c - 1). 

For the 1-loop quantities, Eqs.(47, 48), converting to the representation R is a straight- 
forward substitution: 

N f ^N f .(2t R ) (67) 

and, in addition, for Eqs. (41-46): 

c F — > c R (68) 

Aside from these changes, all algebraic expressions (and the numerical coefficients resulting 
from loop integrations) remain the same. 

A similar reasoning applies to the 2-loop quantities in Eqs. (49, 51, 53): For most dia- 
grams, once their value is expressed as a linear combination of cp, c R N c and cpNf , it suffices 
to apply substitutions (67) and (68). The only exceptions are diagrams containing a gluon 
tadpole [diagrams 7, 14 (Fig.2); diagram 3 (Fig.4); 1-loop diagrams, when expressed in terms 
of a R i', a R p by means of Z gi Za]' In these cases, only one power of c R should be changed to 
c R ; a possible additional power of c F originates from a gluon tadpole and should stay as is. 
This peculiarity implies that, in order to perform the substitutions as described above, one 
must start from the per diagram breakdown of 2-loop results. To avoid presenting a lengthy 
breakdown, we apply, instead, substitutions (67) and (68) indiscriminately on Eqs. (49, 51, 
53); consequently, we must then add a correction term, as follows: 
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ryL,RI f I _ rylj^RV I 

\R~ \c F ^c R , N f ^2Nft R 



4 



9Z 



cr {c r -c f )-[- 4tt 2 ln(</T) - 467.9141661(2) 



(16tt 2 ) 2 "" v "" "" , L "" ~"^ Lf 

+88.7817709(1) c sw + 55.1618942(2) <| w ] (69) 



ryL,RI f I _ rylj^RI 1 I 



+ (1^2)2 Cfi (Ci? ~~ ° f) ' [ " 127r ' ln ( a ^ 2 ) + 708.732752(6) 

+305.48068(1) c sw - 54.495244(2) cg w ] (70) 

r^L 7 RI f I _ rylj^RV I 

^P li? - \c F ^c R , Nf^2Nft R 

+ CR (Ci? ~ Cf) ' [ ~~ 121x2 ln (°^ 2 ) + !089.424358(4) 

-88.7817709(1) c sw + 80.378675(2) t| w ] (71) 

[Actually, the reader could arrive at these results without knowledge of the per diagram 
breakdown, by virtue of the following fact: All 'exceptional' powers of cp cancel out of 
Z^ ,RI , Zg )RI , Zp RI , if these are expressed in terms of the renormalized coupling constant 
clri' . Thus, one may: 

• Express Eqs.(49, 51, 53) in terms of gRi> by means of g Q = (Zg' RI ') g R p , with Zg' RI ' 
in the fundamental representation (Eq.(48)) 

• Apply substitutions (67), (68) throughout 

• If desired, reexpress everything in terms of g Q (using (Z^'^ 7 ') -1 from Eq.(48), with 
Nf — > 2Nft R and cp as is) 

No correction terms are necessary in this procedure.] 
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